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Abstract 

The shift operators for XXX-Heisenberg chain are found. They are formed by local 
Yangian operators and the amplitutes of the eigenfunctions obeying Bethe ansatz for the 
Hamiltonian. The physical implication of the shift operators are also explained. 
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Shift Operators for XXX-Heisenberg Chain 



The XXX-Heisenberg chain(HC) defined by 

- 1 
ff = J^(S j -S J+1 - i ), (1) 

i=i 

is exactly solved by the Bethe ansatz [Q. The model describes nearest- neighbour interacting 
spins, situated on the sites of a periodic lattice N. For J > and J < 0, it corresponds to the 
anti- ferromagnetic and ferromagnetic spin chain, respectively. The simple form of equation (1) 
bellies the rich physical behaviour that it displays, and an understanding of the physics of the HC 
in one-dimension has proved a formidable task for theoretical and methematical physicists over 
the last six decades @ [| g [| § @ § @ (H). Explict form of the eigenfunctions with 
r(r = 0, 1, 2, • • • , N) down-spins (i.e. with the eigenvalue S z = — r) are given by|l[ 

| Vv >= ^ a(m 1 ,m 2 ,- ■ ■ ,m r )cf)(m 1 ,m2,- ■ ■ ,m r ), (2) 

m 1 <m2<---<rn r 

where (j)(mi, m 2 , • ■ • , m r ) represents a spin state with r down-spins on the m^-th (J = 1, 2, • • • , r) 
sites. The coefficients 

r! r ^ 

a(m 1 ,m 2 , - ■ ■ ,m r ) = ^ exp[i(^ 9p k m k + - ^ </»j^,,p n )], (3) 

P=l k=l 

are defined only for the ordering mi < m 2 < ■ ■ ■ < m r , and P is any permutation of the r numbers 
1, 2, • ■ • , r, Pfc the number replaced k under this permutation, and tpkn = —<Pnk- In eq.(j|) for the 
case with r = 

Vo >=|TT •••?>, 

is the vaccum state with all spins up. In the Bethe ansatz eigenfunctions ip r > (r = 1, 2, •, N) are 
always generated from | ^ > by the operator 

[Qr= X] a ( m ii m 2, - ■ ■ ,m r )S mi S m2 - ■ ■ S mr ] | Vo >=| A >, (4) 

mi<m2<-<m r 

that has been well known. On the other hand, from the point of view of the operator method, 
if shift operators of a Hamiltonian (or other observables) known, then the eigen-problem can be 
solved exactly. Inversly, if a Hamiltonian system can be solved exactly, in principle, it must have 
shift operators. Shift operator transforming neighboured level states are discussed extensively in 
the linear systems in quantum mechanics, such as the harmonic oscillator and the hydrogen atom 
JL2"| . It is natural to ask whether such kind of operators (not the cases in eq.(|])) do exist in the 
non-linear models, such as XXX-Heisenberg chain. The purpose of this paper is just to establish 
shift operators transform | Vv-i > to | tp r > in the model. 
We start from introducing the unified shift operators 

JV 

Qt-i = -iJ2 W S^xSk)-, W$ = -W%\ (5) 
j<k 

in | Vv-i > there are already (r — 1) down-spins, i.e. acting —i(Sj X Sk)~ = S~S£. — S^S* on 
ipr-i > will lead to a state with r down-spins: 

Qtr-l I Vv-i >=| Vv >, (6) 
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namely, Q^ r -i transforms the adjacent eigenfunctions specified by r — 1 and r. What we need 
to do is to determine the coefHcents W^} for given states | ipr-i > and | ip r >. It is worth 

mentioning that the operator J = — i Sj x Sk is a Yangian operator. Operator Q^ r _ x is yielded 

j'<* 

by combining W- fe with local Yangian operator Sj x Sjt and summation, hence Q^r-i nas the 
most natural form and is a generalization of the Yangian operator. 

Taking eq.(^) and H \ tp r >= E r \ ifj r > into account, one can verify that 

{[H,Q+ r _ 1 ]=U r ,r-lQt,r-l) I VV-1 >, (?) 

where 

i.e. the energy interval. Therefore in the sense of acting the operators on eigenfunctions of H, 
QJr-i are sm ^ operators of H. 

To find the cxplict form of Q J ^ r _ 1 , we need to determine Wp). We consider the cases 
1, 2, • • • successively. 

(a) r = 1. In this case, eq.(||) becomes 

N 

r(l), 



r 



£o = -*E W H $ x i ^ o >= i ^ > • ( 8 ) 
j<k 



Because 



and 



TV 



ipi >= a(m)<p(m), 



m—l 



[ x S k )~ = SjS z k - S^S* ] |tf>= J(llf> - IU>), 



2 

by comparing the corresponding coefficients of <p{m) of the both sides of eq.(^|), we have 

\\W£1 + ■■■ + W^Vi + W Sm+i + ■■■ + W^ N ] = o(m), (m = 1, 2, • • •). (9) 
Under periodic boundary condition it is well-known: 

N 

a(m) — exp(im#), a(m) = 0, 

m— 1 

with 

2tt ,N s N 

9 = — n; n = 0,±l,---,±(y-l),±y. 

Therefore, one obtain the solutions of eq. (|9]) 

W$ = ^(exp(ij0) - exp(i*e)), 0'. fc = 1, 2, • • ■ , JV). (10) 

There are many ways to transform | ipo > to | ipi >, however as we see later, the form of 
cq. (|l0|) can be extented to more general one transforming | ip r -i > and | ip r >. 
(b) r = 2. In this case, eq.([j|) becomes 

TV 

= -*E x &) _ i ^ >= i ^ >, (11) 
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and 

N 

I 1P2 >= a ( m i^ m 2)4'(mi,m2), (mi,m 2 = 1,2, • • • ,N) 

mi <m,2 

where 

mi m2 

4>(mi,m 2 ) =\ ■■■ I • • • I • • • > 

represents spin-down state that on the m^-th (j = 1, 2) sites. From the known result of the Bethe 
Ansatz, one has 

a(m, k) = Ce mBl e lk62 + C'e m62 e tk01 , (m, k = 1, 2, • • • , N) 

with 

C = exp(z — ), C = cxp(-j— ), 



and 



(p\2 ^1 "2 
2 COt = COt COt , — 7T < 012 < 7T. 

2 2 2 _r 



i.e. C and C depent only on 9\ and 02, not on m\ and 777,2 (or m and k ). 
Noting that 

- x 4)" |lf>= \ |U>; x if !>= |U>, (12) 

and by comparing the corresponding coefficients of <j>{m\ , 777-2) of the both sides of eq.(pl|), we have 
1 1 N 

2 - o(ma)] + - ]T ^flK) + W^?, TO o(m a )] = a{m u m 2 ), (13) 

m/mi ,m2 

with mi, ?Ti2 = 1, 2, • • • , N . 

Now we come to solve eq.(|13[), i.e. to find the solutions of Wmi,m 2 - 
Because 

we have 



a(mi) = exp(imi0), a(ma) = exp(im2^), 



9 -a(mi) = a(mi), i-^-a(m 2 ) = 0, 



i(9 9mi ' 7$ 9mi 

^ -0(1712) = afms), a(mi) = 0, 



or in general 



18 dm\ ' 16 drti\ 

1 d 

-r^^-a(j) = 6 jm a{m), (m,j = 1, 2, • • • , N), (14) 
w am 



where m — x m represents the coordinate of the spin stays on the m— th site of the lattice. 
From eq.rtlj), one can obtain 



1 d 1 d N 



7f)7^ 1 ^ >= TflT^T 51 «M^M = a(m)(P(m), (15) 



7$ 9m 70 dm 

m— 1 

hence, the action of the partial differential operator is clear, when it acts on | tpi >, it will 

pick up the term a{m)<p{m) among | ipx >. 
From the above analysis, if we set 
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ith 

a(k,j) if j > k, 



(2) 

One may argue that why is not a number, but a partial differential operator acting on 



A(j,k) = A(k,j) = 

then eq.([l3j) is satisfied 

,t - 

the coefficients a(J). If Wj^ is a number, then eq.(|l3|) can be viewed as a linear N(N — l)/2 set 

(2) 

of equations. In the procedure for determining the coefficients Wj^, one finds that he must set a 
relation between a(j) and a(j, k),(j, k = 1, 2, • • • , N). However, the relation cannot be valid when 
one substitutes the results of a(J) and a(J, k) obtained from the Bethe ansatz. The fact can easily 
be verified by taking the example in the case with N — A. Therefore, W- fe cannot be numbers, 
but may be operators. 

In the following, we shall derive the shift operator Q21 from physical consideration, so that 
the physical picture for the transformation from ^ > to | ip2 > can be seen clearly. In the former 
state, one of N spins is down, while in the latter, two of N spins are down. The crucial point 
is that, based on | tpi > in which there has already been one down-spin, how we can invert the 
second one so that it can be shift to be | ip2 >• Firstly, we must know clearly on which site the 
spin is down. For a given term 

12 rn N 

a(m)<j)(m) = a (to) |TT • • • I • • • t >, 

it indicates that only on the to— th site the spin-state is down, and its amplitute is a(m). With 
the help of the partial differential operator jg ^ , the wanted term can be picked up from | > 
(see eq.(|l5|)). After this manipulation, we know clearly that on the m— th site the spin has already 
been down, thus the second down-spin can occur on the j— th (j = 1, 2, • • • , m — 1, m + 1, • • • , N) 
sites, respectively Meanwhile, | ^2 > can be rewritten as 

N 



1 ^ >= \ E $ ( 



m-l N 

= E a (-?' m )^0' TO ) + E a(m,k)<p(m,k). (17) 

j—l k— m+1 



Introducing 



m-l N 

T m = \ V TT ,. > + V T~ 

j — l fc— m+1 



since 



T m^u, m ) = J \ ( 2i )( S j x S m)~ }a(m)<f>(m) = a(j, m)(f>(j, m), (j < m), 



T m ^(m,k) = a ^'^ (-2i)(^ m x S k ) ]a(m)0(m) = a(m, k)<f>(m, k), (k > m), 



one obtains 
Define 

we have 



a(m) 

'm, k] 
a(m) 

T m a(m)4>(m) = $(m). 

F — T IA 
w am 

m-l N 

F m \i>i>=2_j + E a ( m ,k)(f>(m,k) 

j=l k=m+l 
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so that 



/V 



(J2 Fm) | V>1 >=2 I V-2 > • (18) 



Since (see eq.(||)) 



we then arrive at 



N 



i w 



m— 1 

by identifying the coefficients of (Sj x Sfc)~ of the both sides of eq.(]l^), it leads to 

r (2) 1 , . x r 1 8 1 d 



WZ = -^a(j,m) 



h ' m i9 ' a(j) dj a(m) dm 

W£\ = ~a(m,fc) [-±-1 (20) 

' ' it) a[m) am a[k) ok 

i.e. they are nothing but the explicit expression of eq.(|l6|). 

Therefore, the physical picture for the transformation from | ipx > to | tp2 > is clear. This idea 
can be extented to generate | ipr > for given | Vv-i >, for example, in the similar manner we can 
find the shift operator transforming \ ip2 > to | ip 3 >. 

(c) r = 3. 

Similar to (b), the eigenfunction 

N 

I V>3 >= ^2 a(m 1 ,m 2 ,m 3 )(l>(m 1 ,m 2 ,m 3 ), (m 1 , m 2 , m 3 = 1, 2, • • • , N) 

mi<mj<iti3 

can be recast to 

1 W 

I ^ 3 >= 3 X! $ (™1> TO 2), 

mi <ma 

where 

m i — 1 



$(mi,m 2 )= a(j,TO 1 ,TO 2 )0(j,TO 1 ,TO 2 )- 



m 2 -l JV 

a(mi,k,m,2)<p(mi,k, 7TI2) + ^ a(m 1 ,m 2 ,l)(t)(m,i,m 2 ,l). 

fe-mi+l /— 7712 + 1 

In I ^2 > two of 2V spins are down, while in | ^3 > three are three down-spins. Just like the 
case in (b), at first, we introduce the partial differential operator j^j^g§^g§^, which yields 

119 5, 

' >= a(mi,m 2 )0(TOi,TO 2 ), (21) 



iOi 162 dmi dm 2 
i.e. the term 

777 1 7772 

a(mi,m 2 )</>(mi,m 2 ) = a(m 1 ,m 2 ) | • • • | ••• | •••> 

is picked up from | ip 2 >, so that two down-spins on the mi— th and m 2 — th sites have been pre-set. 
Further, the third down-spin can occur on the j—th site: 

(J = 1, 2, • • • , mi — 1, mi + 1, mi + 2, • • ■ , m 2 — 1, m 2 + 1, m 2 + 2, • ■ ■ , N). 
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Introducing 

mi— 1 m2-l N 

(mi, 7712) / ; (mi,m2)^(j,mi.m2) / < [rn\ ,m 2 )^(mi ,k,rri2 ) / ^ (mi,m2)^(mi,m2,i) 1 

i=l fc-mi+l I—7-H2 + I 



where 



and 



T7 w . . = a(j ' Wl,m ^ (2i)[ (& x S m J- + (Sj x S m2 )" ]: 

(mi,ro 2 )-t(;,mi,ro 2 ) a{mx,mi) 
T (mi,m 2 )^(mi,fe,m 2 ) = a {^l um A ^ ~^ TOl X + ^ X ^ m2 ^ ^ 

T7 w n = a(mi,W2 \^ (2t)[ -(5 mi x - (SL x ]■ 

(mi,ros)^(mi,m2,l) a(mi UI2 ) 



One can verify that 

T , (mi,m 2 )a(mi,m 2 )0(mi,TO 2 ) = $(mi,TO 2 ). 

further, define 



-F(mi,m 2 ) -T (miim2) — — — — 

we get 



N 

( E ^1,%)) I ^2 >= 3 I ^3 >, 
mi <m 2 

on the other hand, (see eq. (J5]) and eq. 

N 

[ Qt,2 = ^kiSj X ] I V>2 > = | ^3 >, 

thus 

1 " 

^2=3 E F {mum 2 ), (22) 
mi <m,2 

by making comparison the coefficients of (Sj x of the both sides of eq. (|2^) , it yields 

with 

{a>U,k,l) if j <k<l, 
a(l,j,k) if Kj<k, 
a(j, I, k) if j <l < k. 

Thus the operator Q% 2 shifting | -0 2 > to | ^3 > is found. 

Making use of the similar analysis, one can obtain the general form of Q^ r _ 1 with the solution 

of 



.,, _ 2 1 1 



1 9 1 d , d d d 



A(j, h,h, ■ ■ ■ Jr-2) dj A(k,h,h,- ■ ■ Jr-2) dk dh dl 2 dl r - 



}. (24) 
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In (b) and (c) one sees that the actions of the partial differential operators -4 S~ and -4- -4- ?r— 

\ / \ / r 1 iv ova i&i i&2 ami ora^ 

are picking up certain terms a(m)<f>(m) and a(777i, 7772)</>(mi, 7772) from | ipi > and | t/>2 > respec- 
tively. However, since | ipo > has only one single term, we need not to introduce some partial 
differential operators to pick it up from | ipo >. That is the reason why partial differential opera- 
tors do not emerge in the coefficients Wj 1 ^. 

An arbitrary state | ip r > can be obtained by repeated application of Q^ r _ 1 to the state | i/jq >: 

I A >= Q+ r _ 1 Q+_ 1)f ._ a • • • QtiQto I >, (25) 

and a caution should be made that when refer to an operator Q^ r _ 1 , it always acts on | ipr-i >■ 
The facts in eq.(^) and eq.(p5j) are mostly similar to those in a hydrogen atom, where the 
eigenfunctions are expressed in terms of the radial wavefunction R n i(r) and the angular part 
wavefunction Yz m (r) as 

(r) = R n i(r)Y lm (v), (26) 

here r is the three-dimensional coordinate, f = r'r a unit vector, Yj m (r) is a spherical harmonic 
and 

Rni{r)=C n i 1 F 1 (-n + l + l;2l + 2;— )(— )'exp(-— ) 

na na na 

with C n i the normalized constant and a the Bohr radius. 

For a given state V'z+i.u — Ri+i,iYu, the unit dipole vector f will shift the angular part 
wavefunction as follows: 

where f + = x + y. 

To shift the radial part one must introduce the dilatation operator 

£>^ t i, n = ex P [(irp r + -)in-J T ] ) (27) 
h 2 nil 

with the radial momentum p r = —ih(-^- + £), a partial differential operator. Since 

± r _ ^_ D ± 
n±1 ' n n ~ nil n±1 ' T " 

it leads to 

Dt +2>l+1 j^Ri + i,i(r) = -\W + 3)(2I + 4)] 1 /2^ +2ii+l(r) . 

^^iT^Y^^+i.mW = a[(Z + + 2)]VV,+a J i+i,i+i(r), (28) 
and a general state ^j+i,i,j(r) can be generated from the ground state , 0ioo( r ) by 

^ +1 , M (r) = a-\l\)-\l + 1)- 1 / 2 (D+ +U j) ■ ■ ■ (I?^)(I)+^)(r+) Wr). (29) 

Eq.(||) and eq.(|29|) are quite similar to eq.(p8|) and eq.(25), respectively, based on the follow- 
ing sense: for the Bethe state | ip r > (see eq.(|2|)), it also consists of an angular part wavefunc- 
tion <p(mi, 7712 5 ' ' ' i m r) and a radial part (one-dimensional) wavefunction a(mi, m,2, • • • , 777 r ) that 
expresses the phase transition. In the shift operator Q^ r _ 1 , the action of —i(Sj x Sk)~ is trans- 

forming the angular part, and the partial differential operaor Wj k transforms the radial part, just 
like the dilatation operator in eq.((2l) dose. Certainly, there are some differeces between the Hy- 



As a result 



drogen atom and the XXX-Heisenberg chain, for the former its wavefunction ipni m is simultaneous 
eigenstate of the set {H at0 m, L 2 , L z }, while for the latter \ i/j r > is simultaneous eigenstate of the 
set {Hxxx, S z }, but not an eigenstate of the square of the total spins operator S 2 . 
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Eventually, we would like to point out: (i) the shift operator Qi can also be obtained from 
physical consideration the same as that developing in (b) and (c) (see Appendix A); (ii) From 
eq.(^), we see that Q^r-i are raising operators, which can also be found following the same idea 
(see Appendix B). In conclusion, we have established shift operators for XXX-Heisenberg chain. 
It is interesting to extent the method to the Hubbard modelp3 1 114 and the generalized Bethe 
ansatz^] in subsequent investigations. 
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APPENDIX A: Clear Physical Picture for Shifting | ip > to | Vi > 



Because 



jk 1 jk jk jk jk jk 

2< 



■US, x s t y |n>= 5 (|iT> - ITI>), (SJ + s;) ITT>= (HT> + ITi>), (30) 



so that 

3 k j J fc 

|4T>= -i(Sj x S fe )- + -(fir + Sj;) |TT>, 

|U>= x S k )- + \(Sr + 5 fc -) |ff > ■ (31) 
Due to eq.(|3l|), one finds the transformation 

12 m N 12 m N 

ITT ••■ T ••• ]>=^a(m)cj ) (m) = a{m) |TT — I T> (m = 1, 2, • • • , iV), 
can be realized by the operator: 

1 m— 1 iV 

^ m = ^T3T«( m ){ E t x + 2 {S ^ + Sm) ] + E t x + + 5 7) ]>■ 

(32) 

Define 

TV 

^0 = E ^O^m = 
m— 1 



- ^aT^i X>0') - x &)~ + E t^y^M + ^ E a ( fc )l^' ( 33 ) 

j<fe m— 1 fc— 1 

one gets 

N 

F n | >= E a{m)4>{m) =1 > ( 34 ) 

Using 

TV 

E a ( m ) = °' 



and 



iv 



E a ( m ) S m I "00 >=| ipi >, (35) 



from eq. (p3|) and eq. (p4|) we then have 

1 i i\r_o 

" l J^Y EHJ) - <k)](S 3 x S k )~ | Vo >= (1 - g y) I ^ >, (36) 

thus 

2 * 

- z ^ E^j) - °(*)K S i x &)~ I >=| Vi >, (37) 

j<k 
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on the other hand 

JV 

.(1) 



[ Qi = -iY, W * & x ] I V'o >=| V>i >, (38) 

j<k 

by comparing the coefficients of (Sj x S k )~ of both side of eq.(^) and eq. (|37|) , it leads to 

W$ = l(exp(ijO) - exp(ik9)), (j, k = 1, 2, • • ■ , N) (39) 

which is nothing but eq.([Tc|). Thus the physical picture for the transformation from | ipQ > to 
| Vi > is also clear. 



APPENDIX B: The Lowering Operators 



Because i(Sj x S k ) + — Sj~S% — S^Sj, one can calculate that 

^(S f J xS k )+ \U>=-\ |tf>; KSj x S k ) + ||T>= \ |T T> ■ (40) 

We set the Ansatz: 

N 

[Qo,i = iJ2 W h (1) j xS k )+] |Vi>=|V>o>, W,' fc (1) = (41) 
i<fc 

we obtain the equation 

1 AT m-1 AT 

iE[(-E^ (1) + E W; fc (1) Mrn)] = l, (42) 



whose solution are 



Becuase 



2 

m=l j'=l k=m+l 



T , (2) 2 i r i 9 i a 



3 > k N(N — 1) iO a(j) d(m = j) a(k)d{m = k) i ' 1 ' 



i(Sj x S k ) + |I|>= i(|lf > - |U>), (5+ + S+) ||{>= (|lf> + \\t>), (44) 

so that 

|lf>= i0i x 5 fc )+ + + 5+) |U>, 

IU>= x S fe )+ + i(S+ + 5+) |U> . (45) 

Therefore, based on eq.(|45|) there would be some differences in achieving shift operators Q^-i r 
when r > 1. In the position, we set the Ansatz: 

AT AT 

[ or-i,r = *E^ (r) (^ x &) + E A 5 r)5 / ] i ^ >=i ^-i >. ^ (r) = - w U r) - ( 46 ) 
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For instance, when r = 2, since 

rn i m,2 



a(m 1 ,m 2 ) | • • • | • • • | • • • >, 



mi <m 2 



\ _ — , mi m 2 mi m 2 

^ 1>= -^_^ ^ [a(mi)|--- I ••• T •■■ > +a(m 2 ) | •■■ | ■•• I •••> 



mi <m 2 

owing to eq.(Eq), one finds the transformation 



mi m 2 mi m 2 

a(mi,m 2 ) | ••• J. ■•• j •••>=*► a(mi) | •• ■ J, ■•• | •■•> 



can be realized by the operator: 



TTmi^mi = I X Sm 2 ) + + \(S^ + S+ 2 ) ], (47) 



and the transformation 



mi m2 mi 

a(mi,m 2 ) | • •• I ••• I • • • >=^ a(m 2 ) | • • • T ■•• I •> 
is realized by the operator: 

r (mi,m 2 )_ 2 = ^^y[ >< 3n.) + + ^ + ]■ (48) 

Hence, 

[ T (mi,m 2 )^mi + T (mi,m 2 )^m 2 1 «( m l » ^l^K > m 2 ) = «( ffl l ) + «(™2 )<M™2 ) . (49) 

Define 

p' = [T~ + T~ 1 1 1 - - ("501 

(mi,m 2 ) L (mi^H^ ( m ^ m ^ m ^ie 1 i6 2 dm 1 dm 2 1 V j 

so that 

TV 

( E ^(mi,m 2 ))|^2>=(^-l)|^l>, 
mi <m 2 

on the other hand, 



N N 

n ( 2 )/c „c\+V \( 2 ). 
i<fe i 



>r 2 = *E^ (2) (^ x ^) + E A f^ 



thus 

JV 



E ( 51 ) 

mi <m 2 

by making comparison the coefficients of (Sj x Sk) + and S 1 ^ of the both side of eq. (|5l|) , it yields 



w (2) = 1 a(m 1 )-a(TO 2 ) 1 1 d d_ 

mi - m2 N - 1 o(mi,mj) Mi %Q% dm x dm 2 ' (mi m2h 

A (2) = J- 1 1 [ y fl ( m ) ± a ^') 9 9 V a (i) + a ( fe ) 19 9 1 (ko) 
3 N-liQ 1 iQ 2 [ £^ 1 a(m,j) dmdj a(j,k) djdk V 1 ' 

Lowering operators for r > 3 can be obtained in the same way. 
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